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ABSTRACT 

Using a quantitative model for bipolar outflows driven by hydromagnetic protostellar winds, 
we calculate the efficiency of star formation assuming that available gas is either converted 
into stars or ejected in outflows. We estimate the efficiency of a single star formation event 
in a protostellar core, finding 25% — 70% for cores with various possible degrees of flattening. 
The core mass function and the stellar initial mass function have similar slopes, because the 
efficiency is not sensitive to its parameters. We then consider the disruption of gas from a 
dense molecular clump in which a cluster of young stars is being born. In both cases, we 
present analytical formulae for the efficiencies that compare favorably against observations and, 
for clusters, against numerical simulations. We predict efficiencies in the range 30% — 50% 
for the regions that form clusters of low-mass stars. In our model, star formation and gas 
dispersal happen concurrently. We neglect the destructive effects of massive stars: our results 
are therefore upper limits to the efficiency in regions more massive than about 3000 Mq. 



1. INTRODUCTION 



Langer, & Wilson 1986; Umemoto, Hirano, Kameya 



Star formation in molecular clouds is a slow and 



Fukui, Kuno, fc Takakubo 1991; Bally fc Duvert 1994; 



Bally, Reipurth, Lada, &: Billawala 1999]) or from 



seemingly inefficient process. Clouds are transformed 
into stars far more slowl}^ than gravity would allow 



protostellar cores (Myers, Heyer, Snell, & Goldsmith 



( ^uckerman fc Palmer 1(5 3 ); onl}^ the deepest regions 



1988|; pakano, Hasegawa, fc Norman 1995 ; Velusamy 



of clouds participate in star formation (McKee 19811 ^ 



Langer 1998| ; |Ladd, Fuller, fc Deane 1998| ). In 



Li, Evans, fc Lada 1997| ) ; star- forming regions are dis- 
rupted before star formation is complete ( Lada, Mar-| 
gulis, &: Dearborn 1984| ); and individual protostellar 
cores may also shed mass on their way to stardom 
( [Benson fc Myers 19891 ; |Ladd, Fuller, fc Deane l"998|) . 
A growing appreciation for the violence and ubiquity 
of protostellar outflows has led many authors to sug- 
gest that outflows provide the mechanism to regulate 
star formation and to dissipate star-forming gas, es- 
pecially in regions that lack massive stars. Protostel- 
lar outflows have been implicated in the support of 
molecular clouds and cluster-forming clumps against 
gravity (^rman &: Silk 1980|; |Lada fc Gautier 1982| ; 
McKee 1989| ; pertoldi &: McKee 199q ) a nd in the 



removal of gas from star- forming regions (Levreault 
1984 iLanger, Frerking, fc Wilson 1986| ; iGoldsmith 



this paper, we shall construct quantitative models to 
determine the destructive effects of outflows on the 
molecular cores that make individual stars, and on 
the larger molecular clumps that make star clusters.^ 
Disruption by outflows limits the efficiency of low- 
mass star formation; it relates the mass function of 
cores to the stellar initial mass function; and it de- 
termines whether a young stellar system will remain 
bound. 

1.1. Disruption of Protostellar Cores 

The most direct progenitors of stars are dense cores 
within star-forming regions. The distribution of ini- 
tial stellar masses is therefore controlled by two fac- 
tors: the distribution of core masses, and the frac- 
tion of a given core that can successfully become a 



^The nomenclature for the structure of molecular clouds has yet to be established. We adopt the usage of [Williams et al. (200C ): 
Clumps are coherent regions of molecular gas inl — h — v space, generally identified from spectral line maps of molecular emission. 
Star-forming clumps are the massive clumps out of which stellar clusters form. Finally, cores are regions out of which individual 
stars (or multiple systems like binaries) form. We also distinguish protostellar winds, which are ejected from the protostar or the 
protostellar disk, from protostellar outflows, which include material in both the swept-up ambient medium and the wind. 
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star ( |Nakano, Hasegawa, &i Norman 1995| ). Recent 
observations of the core mass distribution by [Motte| 
et al. (19981) indicate a similarity to the initial mass 
function (IMF). 

Several authors have made theoretical predictions 
for the manner in which a protostar's wind truncates 
its own accretion, limiting the core mass that becomes 
stellar. Silk (1995 ) matched the Kelvin-Helmholtz lu- 
minosity of the protostar with the accretion luminos- 
ity of its disk (in the model of Terebey et al. 1984 ) to 
find the point at which disk accretion stops. Sim- 
ilarly, Adams &: Fatuzzo (199^ ) matched the ram 



observations of core disruption and constraints on the 



efficiency of core collapse in §5.4 



pressure of material falling directly onto the stellar 
surface (in the same infall model) with the wind ram 
pressure; although it is more physical to consider ram 
pressures than luminosities, this theory does not ex- 
plain what stops disk accretion. [Nakano et al. (1995 ) 
took a distinctly different approach: recognizing that 
material falling onto the disk might circumvent the 
outfiow, they instead considered the disruption of 
the larger-scale density distribution of the accreting 
core. Idealizing the core material as static, they ar- 
gued that accretion would stop once a momentum- 
conserving wind bubble (or a stellar HII region) had 
traversed its radius. For fiducial values of the param- 
eters, they find quite low efficiencies, of order a few 
percent. 

All of these models have considered only spheri- 
cal protostellar winds, and they have generally not 
treated anisotropics in the protostellar core material 
(besides centrifugal deflection of the infall). However, 
the winds from young stars are collimated into jets 
on a short length scale, and the core is likely to be 
somewhat flattened because of support from poloidal 
magnetic fields. Li Sz Shu (1996| ) incorporated these 
very elements into a model for the early phase of a 
molecular outflow, during which it propagates into 
its own protostellar core. However, Li & Shu did not 
consider the implications of this model for the effi- 
ciency protostellar accretion. 

In §|5| we shall construct a model for the disrup- 
tion of an accreting, somewhat flattened protostel- 
lar core. Assuming that the wind intensity is di- 
rectly proportio nal to the accre tion rate, as in the 
X-wind mode l (^hu et al. 1994 ), and in some disk- 
wind models (pelletier fc Pudritz 1992| ), we identify 
the angle from the axis within which core material 
is ejected and outside of which it is accreted. As in 
the model of Nakano et al. (1995 ), and also as ar- 
gued by Mouschovias (199l| ), protostellar accretion is 
terminated by the flnite size of the protostellar core; 
however, the efficiency with which the star forms is 
a function both of the flattening of the core, and the 
ratio of its sound speed to the characteristic wind ve- 
locity. We compare the results of this calculation to 



1.2. Disruption of Cluster- Forming Clumps 

A protostellar wind can have destructive effects on 
regions considerably larger and more massive than its 
own protostellar core, repercussions that have impli- 
cations for the fates of young stellar clusters. 

Stars are born in groups, deeply embedded within 
dense clumps of molecular gas. The fate of a nascent 
stellar system - whether to become an expanding as- 
sociation or to remain a bound open cluster - is de- 
termined by how rapidly and effectively it breaks free 
of the gas from which it was created. This process is 
primarily controlled by two parameters: the efficiency 
of star formation, defined as the fraction of gas con- 
verted into stars, and the time scale on which the re- 
maining gas is disrupted. If the gas removal is rapid 
relative to the free-fall time tg, then more than half 
the mass must be in stars for the cluster to remain 



bound ( Hills 1980 ). On the other hand, sufficiently 
slow mass loss allows a virialized stellar system to 
expand adiabatically and to remain bound. [Lada 
|et al. (1984| ) studied intermediate cases, in which they 
found that part of the cluster was likely to escape even 
if a fraction remained bound. 

Empirically, roughly only a tenth of all stars re- 
main bound in clusters ( Miller fc Scalo 1978 ). The 
interactions between stars and gas must produce this 
statistic, yet the underlying mechanism by which gas 
is dispersed remains largely mysterious. Most stars 
are born in large OB associations (e.g., McKee & 
Williams 1997), and the violent effects of the mas- 
sive stars in these associations have the potential to 
disrupt the star- forming gas; however, this can only 
be in addition to the disruptive effects of the low-mass 
stars themselves. 

We shall take the embedded clusters in L1640 as 
prototypical examples of groups of low-mass stars in 



the process of formation. As described by Lada et al 
1991| ) and Lada (199^ ), these clusters have several 



hundred stars apiece; they are coincident with the 
densest and most massive molecular clumps in the 
region. Stars currently constitute about 40% of the 
total mass in the two largest clumps, but the oth- 
ers show a much smaller stellar fraction. Moreover, 



star formation is restricted to these regions (Li et al 



1997|). This obse rvation is consistent with the theory 
of pVIcKee (1989 ) that stars form only in regions of 



high extinction and therefore low ionization, so that 
ambipolar diffusion can lead to the creation of proto- 
stellar cores. It lends credence to the idea that the gas 
and stars can be treated finite, self-gravitating 



system of a definite mass and radius (see also Mathieu 
1983), a notion we shall use in our analysis. 
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In §|6|, we calculate the amount of gas removed from 
a clump by the bipolar molecular outflows driven by 
the stars themselves. There is ample evidence for 
the importance of gas removal by outflows from star- 
forming regions: note, for instance, the disruptions 
of PV Cephei (|Levreault 1984|), B335 ([Langer et al 



1986|)|, Lynds 1221 (pJmemoto et al. 199 1|), HH83 



( iBally fc Duvert 1994| ), and Circinus ( |Bally et al 
1999|]|, and the eruptions of the IRS3 outflow from 
Barnard 5 ([Goldsmith et al. 1986|) and of the H H300 
outflow from Taurus ( [Arce Sz Goodman 2000| ). In- 
deed, any example of protostellar jets and Herbig- 
Haro objects outside the source's parent cloud implies 
that potentially star-forming gas has been ejected. 
Individual protostellar winds are thus remarkable in 
their ability to affect pre-stellar cores (on ~ 0.05 pc 
scales), cluster forming clumps (~ 0.4 — 1 pc) and 
giant molecular clouds (~ 10 — 100 pc), as evinced 
by the existence of parsec-scale flows ( [Reipurth et al. 
1997|]. 

We shall gauge mass ejection using the model for 
protostellar outflows presented in a previous paper 
( [Matzner fc McKee 1999| ). This model generalizes 
the force distribution predicted in the X-wind model 
( ghu et al. 19951) to a much wider class of winds. 
The interaction of the protostellar wind with the am- 
bient medium occurs in a thin, radiative shell that 
conserves the wind momentum and ambient mass in 
each direction. In §Q we use this model to calculate 
the criterion for the outflow shell to be ejected in a 
particular direction. The mass loss from such out- 
flows is estimated analytically in §p. The principal 
idealization in our treatment is that each star forms 
at the exact center of the clump distribution. There- 



fore, we present in §3.1 a numerical evaluation of the 
gas lost from a model clump, in which star forma- 
tion is distributed according to the local ambipolar 
diffusion rate ( McKee 1989| ). In the Appendix, we 
evaluate the detailed effect of gravity on the motion 
of an outflow shell. 

Our analysis suggests that a significant amount of 
mass can be lost during the process of star forma- 
tion itself. Indeed, the rate at which mass is lost due 
to outflows is several times the rate at which gas is 
converted into stars. However, we do not consider 
the disruptive effects of massive stars, which might 
cause gas to be removed more rapidly. Our results 
are therefore only rough upper limits for the star for- 
mation efficiency in clumps more massive than about 
3000 M0, in which O stars will form. 

Unlike most previous work (e.g., Nakano et al. 
1995|]|, we divide the calculation of star formation ef- 
ficiency into two parts: the disruption of a star's own 

^We use lowercase (m*, m^) for the mass of an individual 
cloud mass or the mass of an ensemble of stars. 



protostellar core, and the disruption of a clump in 
which an aggregate of young stars is forming. Dif- 
ferent assumptions are appropriate for the two cases. 
For a core, the rates of accretion and outflow are di- 
rectly proportional; the axes along which the core 
is flattened and the wind is collimated are likely to 
be aligned; the wind duration exceeds the free-fall 
time; and the density distribution declines with ra- 
dius approximately as r~^. For the disruption of a 
cluster-forming clump, opposite assumptions apply: 
for instance, winds are impulsive and not likely to 
correlate in direction with the clump's structure. 

1.3. Estimate of the Star Formation Efficiency 

Before getting into the details of the calculation, it 
is worthwhile to make a simple estimate of the out- 
come. The star formation efficiency for the formation 
of a single star is 



(1) 



where Mej is the mass^ ejected from the protostellar 
core or star forming clump due to the formation of 
a single star of mass (see §^ below). The total 
momentum in the wind is = rriwVw. We assume 
that the mass of the wind is a fraction fu, of the stel- 
lar mass, so that = fwf^* and Pw = fwfn^Vw 
This momentum is transferred to the cloud by ra- 
diative shocks; as described above, such a model is in 



good agreement with observation (Matzner & McKee 
1999| ). In order for the shocked gas to escape from 



the cloud, it must be moving at least at the escape 
velocity, fesc- Since the shocked gas has momentum 
Pw, the ejected mass is a maximum if all the gas has 
this minimum velocity, Mej ~ Pw/vesc- Assuming for 
the moment that this mass is large compared to the 
mass of the star that forms, as in the model of |Nakano| 
et al. (19"95| ), we estimate the star formation efficiency 



as 



Pv 



fwVu 



(2) 



where the factor Cg > 1 allows for the decrease in the 
momentum of the shocked shell due to the self-gravity 
of the core or clump. 

This argument can be refined by taking into ac- 
count the central fact about protostellar winds, that 



they are collimated. As originally shown by Shu et al 
1995|) and demo nstrated more generally by Matzner 



McKee (1999| ), the wind momentum falls off ap- 
proximately as l/(sin0)^ for 9 > where 6 is mea- 
sured from the axis of the wind and 9o ~ 10~^ is 
discussed in §^ below. The momentum inside 6 thus 
scales as ln(l/0): equal amounts of momentum are 
star and its wind, and uppercase (Moj, M^i and M*) for ambient 
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contained in each octave of angle. Since the area cov- 
ered by the wind, and therefore the mass swept up, is 
dominated by the largest angle, it follows that most 
of the mass that escapes will indeed have the min- 
imum possible velocity (~ Vesc), as assumed above. 
Furthermore, since the momentum is spread out log- 
arithmically, only a fraction ~ l/ln(l/0o) of the total 
momentum is available to eject matter at ~ Vesc- As 
a result, the star formation efficiency estimated above 
should be multiplied by this quantity. The calcula- 
tions below show that this fraction is more precisely 
given by ~ l/ln(2/0o)- Defining the efficiency fac- 
tor X by 

'"P/""!, (3) 



we conclude that for small values of X, the star for- 
mation efficiency e ^ X] sX large values of X, we 
expect e to approach its maximum possible value, 
unity. These expectations are borne out by the re- 
sults presented below. Observe that the collimation 
of the wind reduces the mass that can be ejected by 
a factor ~ l/ln(2/0o) ~ 0.2 and increases the star 
formation efficiency by the same factor. As a result, 
our estimates of e are significantly higher than those 



of Nakano et al. (1995 ) and are in much closer agree- 
ment with observation. 

2. INSTANTANEOUS, OBSERVATIONAL AND TOTAL 
STAR FORMATION EFFICIENCIES 

We begin by discussing the various possible defini- 
tions of the star formation efficiency. For the collapse 
of a protostellar core to form a single star (or a bi- 
nary) there is only one efficiency that can be defined, 
namely the ratio of the star's (or binary's) mass to 
that of its parent core: 



(4) 



This definition is consistent with that in equation 
because when the star formation process is complete, 
all the mass has been incorporated into the star(s) or 
ejected from the system. 

For the formation of a stellar cluster in a molecular 
clump, we generalize e so that it is the instantaneous 
star formation efficiency, which is the rate at which 
mass is converted into stars, divided by the rate at 
which mass is lost from the gas. (Implicit in this def- 
inition is the notion that the star-forming clump is a 
well-defined gaseous system.) We assume that only 
two processes change the gas mass: star formation, 
and mass ejection due to outfiows. Then, the instan- 
taneous star formation efficiency is 



dMci dM^ + dMe 



(5) 



where is the mass of a set of stars, is the 
mass of the clump, and Mgj is the total mass ejected 
from the clump by protostellar winds; we anticipate 
that part, if not all, of the mass in the winds will be 
ejected from the clump, and this mass is included in 
Mej . For the formation of a single star, dm-), becomes 
and dMej becomes Mef, result, in this case 
e is the same as that in equation (jl]) and is identical 
with Ecore- Becausc of this generality, we shall often 
refer to M^i as the "cloud" mass, using "cloud" to re- 
fer to either a clump or a core. As we shall see below, 
in our model the ejected mass is directly proportional 
to the stellar mass. As a result, e can be evaluated 
at any convenient mass, such as 1 Mq: 







M0 + Mej(l) 



(6) 



where Mej(l) is the mass ejected from the cloud due 
to the formation of a star of 1 Mq . If gas is disrupted 
by other means, as might occur where massive stars 
are forming, then equation (^) will be an overestimate 
of e. 

The observational measure of star formation effi- 
ciency, SFE, is the current fraction of the total mass 
that is stellar: 



SFE 



+ Mel 



(7) 



This quantity varies from zero to unity as a group 
of stars forms and is revealed. Under the assump- 
tion that all the stars that have formed in the clump 
can be identified, even if the stars are unbound, the 
mass of the clump evolves as M^i = Md^Q — J e~^(i-M^, 
where M^i^ o is the clump mass before any stars form. 
As a result, we have 



SFE= i + ^§^--^rsW 

I MWo 



\ /Ah 



-1 



(8) 



The total efficiency of the star formation process, 
SFEtot , is the fraction of the initial clump mass M^^ q 
that will ever become stars. Star formation is fin- 
ished once there is no more gas, and SFE = 1. From 
equation (^), we see that this occurs once (e~^)j^/ = 
Md,o/M*; therefore 



SFE 



tot 



/Ah 



(9) 



where the average is taken over the entire star forma- 
tion process. The instantaneous and total efficiencies 
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are approximately equal, so long as e is relatively 
constant. 

On the other hand, e can be considered a typical 
value for the observed SFE only if both are small. 
For instance, consider a clump in which e is con- 
stant and star formation is half complete (M^ = 
SFEtotMci,o/2 = eMci,o/2). Equation <@ shows that 
SFE = 1/(1 +£-1) at 'this point. 

3. STRUCTURE, INTENSITY AND DURATION OF 
PROTOSTELLAR WINDS 

The angular distribution of the momentum in 
coUimated protostellar winds can be written as 



dpw _ Pwjt) 



(10) 



where /i = cos 9 labels directions relative to the wind 
axis, and P{n) is a force distribution normalized such 
that fn PifJ-)dii = 1. [Matzner fc McKee (19"99| ) and 



Matzner (1999| ) showed that any hydromagnetic pvo- 



tostellar wind should limit to a common force distri- 



biitinli on sralps largp rnmparpd to thp sniirnp- 



1 



ln(2/^o)(l + ^o'-^') 



(11) 



with ^0 ^ 1- Moreover, under the assumption of a lo- 
cally momentum conserving interaction with the am- 
bient cloud ( Shu et al. 1991| ), this relation was shown 
to give good agreement with the observed position- 
velocity and mass-velocity distributions of bipolar 
molecular outflows. In cases where we are interested 
only in angles away from the axis {1 — fi'^ ^ 6q), such 
as the disruption of a protostellar core by its wind 
©), equation (11) simplifies to, 



1 



ln(2/6'o)(sinl 



(12) 



The dependence pw oc (sin 6) ^ is caused by the bal- 
ance of magnetic stresses in the wind at large dis- 



the protostar that drives the wind. Therefore the net 
wind momentum is 



Pv 



(13) 



We shall assume for simplicity that fwVw does not 
depend on m^. This assumption is justified in the 
X-wind model ( |Shu et al. 1994 ) by the insensitivity 
of both and to the details of the accretion. In 
this model, is set by the geometry of field lines 
in a very small region around the disk truncation ra- 
dius. Also, deuterium burning introduces a thermo- 
stat that maintains the escape velocity from the stel- 



lar surface near a constant value ([Stahler 1988]) , and 
the ratio between v^j and the escape velocity is insen- 
sitive to the star's mass, magnetization and accretion 
rate. [Matzner et al. (2000| ) argue that the model of 
Najita fc Shu (1994| ) implies f^Vu, ~ 40 km s"^ if 
one takes the average wind luminosity to be half its 
final value. For a wind velocity of 200 km s~^, this 
implies fw ~ 1/5, between the prediction s ~ 1/3 
and fyi] ~ 1/10 of the ori ginal X-wind (Shu et al 



19881) and disk wind (e.g., pelletier fc Pudritz 1992"|) 



theories, respectively. 

The wind momentum output may vary rapidly, but 
it declines in intensity over a wind or accretion time 
scale tw, which is generally taken to be about 10^ 
years. There is little observational evidence that 
varies with m^,; Bontemps et al. (1996| ) argue for an 



exponential decline in outflow intensity, on a time 
scale of 0.9 x 10^ yr, throughout their sample of Class 
and I protostars. Thus, for star-forming clumps less 
dense than 10^ cm~^, whose free-fall times are longer 
than 10^ years, it is clearly a good approximation to 
consider the deposition of an outflow's momentum as 
impulsive {tw < %). 

What happens at higher densities? The duration 
of an intense protostellar wind is identical to the time 
scale of accretion onto the protostar. This in turn is 
essentially identical to the duration of accretion onto 
the combined star-disk system, for a stable accretion 
disk cannot hold most of the total mass ( [Adams et al 



tance rTrom its source, as in the model of phu et alj 1989 ; 3hu et al. 1990[ ). Finally, this accretion time in 



(1995 ). Deviations from this power-law relationship 
can be caused by variations of the conserved vari- 
ables across streamlines (e.g., because of differences 
in initial disk radius), but these deviations are typi- 
cally quite small. The broadening angle could be 
the result of a wandering of the jet, internal shocks 
from variations in the wind velocity, or hydromag- 
netic instabilities; Matzner fc McKee (1999|) esti- 



a magnetized core is roughly twice the thermal cross- 
ing time of the protostellar core if ambipolar diffu- 
sion is not considered 

thus, the accretion 



[ ([Li &: Shu 1997] ), and less if 

it is ( [Safier et al.T997| ; p:.i 1998| ); 



mated Oq ~ 10~ on the basis of comparisons with 
observations. 

The wind momentum p^j is the product of the wind 
velocity v^j ~ 200 km and the wind mass m^, 
which itself is a fraction of the final mass of 



and dynamical times in a protostellar core are simi- 
lar. But, the dynamical time of a protostellar core is 
shorter than that of a molecular clump if the former 
is an overdense region within the latter. Therefore, 
we shall assume that protostellar winds are always 
impulsive events compared to the dynamical time of 
the parent star-forming clump, even for mean clump 
densities above 10^ cm~^. This assumption, which re- 
mains to be tested by observations, will simplify our 
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analysis of the efficiency of stellar cluster formation 



m 



4. MASS EJECTION BY PROTOSTELLAR OUTFLOWS 

To estimate the amount of mass ejected from ei- 
ther a protostellar core or a star-forming clump we 
approximate the shocked wind and cloud material as 
a thin shell whose motion is purely radial. In other 
words, neither mass nor momentum is mixed between 
angular sectors, and there is no lateral momentum 
generated in the shell. These approximations were 
introduced by phu et al. (1991 ), and were shown by 
Matzner fc McKee (1999 ) to reproduce the observa- 
tional features of molecular outflows. For the case 
in which the ambient gas is not isotropic, we shall 
consider only the monopole component of gravity; 
fc Shu (1997 ) demonstrated that this is a viable ap- 
proximation, even in the extreme case of a disk-like 
core. 

Our approach is as follows. A fixed fraction 
of the accreted material is ejected at velocity in 
a wind with angular dependence given by equation 
(PI). Assuming thin, radiative shocks, purely ra- 
dial motion and monopole gravity, we shall identify 
the direction (characterized by the direction cosine 
//esc = cos^esc) in which the wind is marginally able 
to drive the ambient material outward at the escape 
velocity of the cloud. Only material further toward 
the equator can be accreted, since the rest is ejected. 
And, only this accreted material is available to be 
reprocessed in the wind. The fraction of core mate- 
rial that successfully accretes will thus be determined 
self-consistently. 

To simplify the discussion of the condition for a 
protostellar wind to eject gas from the cloud, we shall 
defer including the effect of gravity on the propaga- 
tion of the shock in the cloud to the Appendix. In 
the absence of the perturbation introduced by grav- 
ity, the shell's velocity Vg is set by the total emitted 
wind momentum and the ambient mass swept up in 
each direction in the direction specified by /x: 



dpwjt) 

dn ' 



(14) 



Here we have assumed that the mass and flight time 
of the wind can be neglected in comparison with the 
ambient mass and age of the outflow. This is a good 
assumption in directions in which Vg ~ feso since the 
mass of the wind must be negligible if it is to drive a 
shock at a velocity small compared to the wind veloc- 



ity (Matzner & McKee 1999). We have also neglected 
the effects of the ambient pressure or magnetic field, 
which should be a reasonably good approximation so 
long as the shock is traveling faster than the local 
signal speed. Note that these perturbations are sim- 



ilar in magnitude to gravity if the clump is in virial 
equilibrium. 

We express the mass distribution of the ambient 
gas as 

dMci 



1 



Q{fi)M,i{r) 



(15) 



dn 4tt' 

where Q is normalized by Jq Q{^)d^ = 1 and M(r) 
is the total mass inside r. Since dpw/dn oc P{^i) (eq. 
|T0|), the equation of motion for the shell becomes 



Q(/i)Mciw. = P{n)p^{t). 



(16) 



The condition for a sector of the shocked gas to es- 
cape from the gravitational field of the cloud in which 
the protostar is embedded can be estimated by set- 
ting pw{t) equal to the total momentum of the wind, 
mwVw = rni^fwVw, and by setting the shock velocity 
equal to the escape velocity, Ugsc- The effects of grav- 
ity on the propagation of the shock are discussed in 
the Appendix; they can be incorporated by including 



a factor Cg > 1 that is evaluated there: 



PjfJ-csc) 

CSC J 



m^fwVw = CgMclVesc- 



(17) 



For impulsive winds (those that last a time that is 
short compared to the free-fall time of the ambient 
cloud, tw "C tg), the factor Cg is close to unity. For 
steady winds (t^ ^ %), Cg is proportional to the ra- 
tio t^/tfi: only the momentum transferred during a 
free-fall time is effective in ejecting cloud material. 

In deriving this escape condition, we have assumed 
that the cloud material does not have a significant 
infall velocity in the region of interest {6 < ^esc)- In 
the theory of ^hu (1977 ), the core material is com- 
pletely at rest outside an expansion wave moving at 
the thermal velocity, a. Since fcsc > «, the outflow 
shell encounters no infall whatsoever in this direction. 



However, more complete theories (Safier et al. 1997; 
Li 199^ ) and recent observations (Lee et al. 1999) sug- 
gest an inflow velocity intermediate between the am- 
bipolar drift velocity and the isothermal sound speed. 
This is still significantly slower than Vesc, justifying 
our approximation that material upstream of the out- 
flow shell is motionless in the direction of interest. 

To evaluate the escape condition (eq. |^) , we intro- 
duce the efhciency parameter X defined in equation 
(^) and use the expression for P{fJ-) in equation (|TT| 



{1 + 9q- /iesc^)Q(Mesc) 



(18) 



Numerically, we have 



X = 0.132c. 



ln(2/0o) 



In 200 V 1 km s 



40 km s 

fw^w 



(19) 
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The mass ejected from the cloud by the protostel- 
lar wind is the sum of the swept-up cloud mass inside 
^esc and the wind mass inside ^esc^ 



Mej = Mci / + (1 - </)^)m^, (20) 



where 



p{^Ji)d^i 



(21) 



is the fraction of the wind momentum retained in the 
cloud; we assume that Vw is independent of and also 
that the shocked wind is dynamically cold, in which 
case (t)w is also the fraction of the wind mass that is 
retained. 

The value of the ejected mass determines the in- 
stantaneous star formation efficiency (eq. 0), 



1 Mei 

- = 1 + _SL 

e rUi, 



1 + 



cl 



m 



(22) 



Q{^^)d,j, + (l-^^)f^. (23) 



Inserting equation (Jlq), we find 



1 1 

e = ^ + X(l + ei - /.e.c2)g(Mesc) ^^^^"^^ 



+ (1 - (l)w)fu 



(24) 



If the cloud is isotropic [QifJ-] = 1], the ejected 
mass is simply Mej = Mci(l - /iesc) + (1 - (l)w)mw 
Neither cloud nor wind mass is ejected if //esc = Ij 
since = 1 if the outflow does not reach the surface 
of the cloud; according to equation (|l8|), this occurs 
at a stellar mass 



(25) 



For X < 1, 6*0 ;S 10"^ and Md < 10^ M©, it follows 
that m^^o .^^ 0.1 Mq, so that most stars are more mas- 
sive than rrii^fi. Equation ( |l8|) implies that the mass 
ejected from an isotropic cloud is 



X{1 + fi^ 



+ (1 - </'«,)/«,m^ (26) 



for m^, > m^fi. Note that if ^ rri-kfi, then the 
amount of mass ejected from the cloud is directly 
proportional to the amount of mass converted into 
stars. The resulting star formation efficiency for an 
isotropic cloud is 



1 



- ~ IH 

e X{1 + /iesc) 



1 



) + (1 - 0^)/^ (27) 



for > m^^O) for ^ "i*,0) we have e = 1. Weak- 
ening the wind or increasing the escape velocity (i.e., 



increasing X) increases the stellar mass m^^o that 
cannot eject any material. For a star with > m^,^Q, 
it causes less material to be ejected: hence iiesc in- 
creases with X (as in equation |^^). This also re- 
duces the amount of wind mass lost: therefore, (/>u, 
increases with X. As a result, we see that the effi- 
ciency e is an increasing function of X. 

5. CORE DISRUPTION AND THE EFFICIENCY OF 
SINGLE STAR FORMATION 

We now consider the disruption of a protostellar 
core by a single protostar. The results of the previous 
section remain valid, but we now have an additional 
condition: all the mass in the cloud with > ^esc ei- 
ther goes into the star or into the wind that escapes: 



+ m^(l 



Mr. 



Qin)dfi. (28) 



For a protostellar core, ^esc is not small, as we shall 
see below. Since the wind is strongly concentrated 
toward the axis according to equation (|TT|), it follows 
that the fraction of the wind material that is retained 
by the core is generally small, (/)^ ^ 1. We therefore 
neglect 4>uj and find the star formation efficiency of 
the core to be (ffcore = s): 



1 



(l + /« 



Q{fi)d^i. 



(29) 



5.1. Isotropic Protostellar Cores 



For the case of an isotropic core, equation (p^) re- 
duces to m^(l + /^) 



M, 



core/^esc) 



where we have 
set (pw = as discussed just above. The mass of 
a protostar in a core is always much greater than 
rrii^fi: By definition, could be of order m^^o only 
if /iesc — 1; but then ~ Meore since /t„ < 1, so 
that m^/mj,^o — ^/X9q ^ 1, a contradiction. Since 
Score = we can use equations (p7|) (neglecting m^^o) 
and (p9|) to obtain an equation for Hesc- 



1 + 



1 



A'esc X{1 + /«,)(! + /ie 

solving for /iesc and Score, 



(30) 



- 1 + - 1 + [1 + 4(1 + /.)2X2]l/2 ^^'^ 

Thus, the efficiency of star formation in an isotropic 
core that forms a single star is always less than the ef- 
ficiency parameter, Seore < X; moreover, it is always 
less than 1/(1 + /^) since the wind mass is easily 
ejected. 



2X 
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5.2. Self-Similar Magnetized, Isothermal Cores 



Equation (|31|) provides a general expression for the 
efficiency of star formation in an isotropic protostel- 
lar core. It depends on the properties of the core 
through the parameter X oc CgVesc- To determine nu- 
merical values for the star formation efficiency, and 
to generalize to non- isotropic cores, we must adopt a 
specific model for the core. A class of models that is 
particularly suitable for this purpose is the class of 
magnetized, isothermal core distributions considered 
by Li &: Shu (1996| ). In these models, magnetic sup- 
port causes the density distribution to be flattened 
relative to the singular isothermal sphere, and it in- 
creases the mean density by a factor (1 -|- Hq): 



pir,f,) = il + Ho)Q{fi) 



2^Gr2 ' 



(32) 



where a is the isothermal sound speed (or effective 
sound speed in the case of turbulent support). ^ 
Shu (1996) have calculated Q{fJ-) and (1 -|- Hq) for 
various values of the dimensionless mass-to-flux ratio. 
The results of these calculations can be approximated 

(33) 



where 



kn = 9AH, 



5/4 
• 



(34) 



This fit was made for the values of Hq presented by 
Li &: Shu (1996| ), and is therefore strictly valid only 
for Hq < 1. It was chosen to make the approxima- 



tion of equation (33) fit the correct Q{9) best near 
the equator for each value of Hq, because we expect 
the boundary between accreted and ejected material 
to be near the equator. This fit is accurate to within 
15% wherever it predicts Q{0) > 0.25; for all values 
< Hq < I, this condition is satisfied for 9 > 45°. 
Close to the equator, this fit is good to within several 
percent - except that when Hq = 1, it is 8% below 



the Li & Shu value 



The mass enclosed within a radius r of a Li-Shu 



core is 



M^orc{r)=2{l+HQ)—r. 



(35) 



Li fc Shu (1997| ) have argued that the mass accre- 
tion rate in the self-similar collapse of such a core is 
always within 5% of the value 



ma 



(1 + ^0 



G 



Comparing the accretion rate to the mass in the core, 
we see that the star-disk system swallows all of the 



mass initially within a radius r = at/2 at a time t af- 
ter the onset of accretion. Equation (3(:) should prop- 
erly be considered a lower limit to the accretion rate, 
as it assumes perfect coupling between mass and flux; 
ambipolar diffusion would allow for faster accretion 
( ^afier et al. 199^ ; Li 199 j ). However, as long as the 
core is not highly magnetized {Hq ^ 1, as expected 
from observations; p^i &: Shu 1997|), the idealization 
of strong coupling is a good one ( |Li 1998| ) . 

Equation (36) assumes that mass accretion occurs 
in all directions. However, material within an angle 
^esc of the axis is ejected rather than accreted. The 
gas outside fiesc is the sum of the star mass and the 
wind mass (eq. ^), so that 

JgJ^^ {sin 9)^0+^ d9 
jQ^\sm9)>^o+id9' 



^core(l ~l" fw) 



Q{lJ.)dfi 



(37) 

(Note that we have set (j)^ = m this equation.) Ac- 
counting for this division, the mass accretion rate is 
somewhat smaller than cited above: 

^3 



,{1 + U){1 + Hq)-. 



(38) 



This expression is based on the simplifying assump- 
tion that the wind does not slow accretion in any 
direction that it does not eject material. Since the 
wind mass is related to the star mass by = fwm-i, 
and since the two share the total accreted mass 
('TT-acc = niw + m^), we have = macc/(l + fw) 
and rhw = macc/i«/ {fw + 1)- 

To determine the star formation efficiency, we must 
evaluate the efficiency parameter X for the Li-Shu 
cores. We first determine the factor Cg that incorpo- 
rates the effects of gravity on the propagation of the 
wind-driven shock. The wind is impulsive or steady 
depending on the value of t^/tg. As remarked above, 
mass inside a radius r is accreted at a time 2r/a; un- 
der the assumption that the core is terminated at a 
radius -Rcorei the wind lifetime is then 



2Rr 



(39) 



Rr 



(40) 



By comparison, the free-fall time in the monopole ap 
proximation is 

_ / 37r y /2 ^ vr 
~ [32G^) ~ 4(l + i7o)V2 

Since tw is several times tg, the wind may be consid- 
ered steady. We make the conservative assumption 
that the wind decouples from the shocked shell once 
it leaves the core, so that Cg is given by equation 

(i.e.. 



(36) ( |A19D . Because the Li-Shu core has p 



oc r 



-2 



kn 



2), this becomes 



1/2 



(41) 
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Noting that Wesc/o = 2(1 + Hq)^^"^, we obtain 

a 



X = 12(l + Fo)ln(2/0o) 
= 0.32(1 +/7o) 



ln(200) V0.2kms-i 
40 km s"^\ 

fw'^w I 



(42) 



The star formation efficiency (eq. |37| ) and the es- 
cape condition (eq. |l^) together determine the escape 
angle 9esc- 



{sin 6)^0+^ de 
'(sin0esc)^«+2 



(1 + U)X. 



(43) 



We find that the solution to this equation gives a star 
formation efficiency of the form in equation ( |3l|) with 
the replacement X ^ {1 + Hq)^/'^X: 



2{l + Hof/'^X 



l+[l+4(l + Fo)3(l + /J2x2]l/2' 



(44) 



The vahdity of this approximation, depicted in Fig - 
P implies that an increase in mean density by a 



ure _ 

factor of (I + Hq) is equivalent, as far as the efficiency 
is concerned, to a weakening of the momentum out- 
put in the wind by the factor (1 + Hq)^^"^. One of 
these factors of (1 + Hq) is contained in X; it arises 
from the increase in escape velocity compared to a 
fixed value of v^, as the mean density is increased. 
The other one and a half powers of (1 + Hq) arise be- 
cause of the increasing concentration of infalling ma- 
terial toward the equator, in the direction the wind 
is weakest. In the limit of a completely disk- like core 
[(1 + Hq) » 1], outflow occurs in all directions but 
along the equator (^esc — ^/2), but only the wind 
mass is ejected [^core ^ 1/(1 + /«;)]• 

5.3. The Initial Mass Function 

Having determined ecorei it is possible to relate the 
IMF to the mass distribution of protostellar cores 
( [Nakano et al. 1995| ). Let dA/"* be the number of 
stars formed with masses in the mass range m^, to 
-|- dm-),. Using a similar notation for the proto- 
stellar cores, we assume that they have a power-law 
distribution: 



dAfc, 



din Mr, 



Let 



din Er 



dlnMc, 



(45) 



(46) 



If we approximate as constant, tlien £q 
m^/Mcore OC il^^e- The IMF is then 

dJ\f* cWcore 



d In nii, d In rri-, 



oc m- 



(47) 



To evaluate e', we assume that the flattening pa- 
rameter Hq is independent of the core mass. Then 
Score depends on X, and it is straightforward to show 
from equation ( |4^ ) that dlnecoro/c^liiX < 1, and ad- 
ditionally that dlnecoro/c^liiX is a decreasing func- 
tion of both X and (I + Hq). From the definition (eq. 
1^, we see that X depends directly on the properties 
of the core only through Vesc- In principle, there could 
be an implicit dependence through f^Vw, but we have 
argued in §^ that this is a constant for low-mass stars. 
Since the protostellar cores are gravitationally bound, 
their internal velocity dispersion is a fraction of their 
escape velocity, a oc v^sc- Here, a includes both ther- 
mal and turbulent motions. Since protostellar cores 
are never colder than 10 K, a > 0.2 km s~^: clumps 
must lie to the right of the circles in figure § There- 
fore, Score cannot be less than 0.3 and d\necore/d\nX 
cannot exceed 0.75, the values for an unflattened core 
with 0" = 0.2 km s~^. Now, protostellar cores o bey a 

( Lar- 
)-^/^, we have 



line-width size relation of the form a oc -R^ore 



son T98TI ). Since V^sc oc {Mcore/R c 



X oc Vesc OC Mclr'e^~^'^\ SO that 



din En 



dinX 



dinX dinMc, 



< 0.75- 



dinX 

d in Mrnr 



0.75- 



p 



2p + l 
(48) 

If the motions are primarily thermal, as we assumed 
in ^5.2| , then p ~ and therefore e' ~ 0. On the other 
hand, if the motions are primarily non-thermal, then 
p ~ 1/2 dCaselh fc Myers 1995| ) and e' < 0.19. As a 
result, for all values of Score; and regardless of whether 
the line width is primarily thermal or nonthermal, we 
have 

< e' < 0.19. (49) 

so that dAf^/dinrUi, oc m^, ^'^'^^ i)°<;orc^ Since in prac- 
tice the thermal contribution to the line width is sig- 
nificant for low-mass cores, and since din Score /(^ In X 
decreases as a (and thus X) increases, which reduces 
s' for higher mass cores, we expect s' ^0.1 in general. 
Then, 



am 

Thus, over the mass range in which the core mass 
function is a power law and the star formation ef- 
ficiency is determined by outflows, the slope of the 
IMF is expected quite generally to be just slightly 
flatter than than that of the protostellar cores. This 
theory does not address the IMF of high-mass stars, 
nor of cores sufficiently turbulent for fragmentation 
to occur during their collapse. 
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Fig. 1. — Core disruption and the efficiency of single star formation. Top: the angle Ocbc (equation [liSj) dividing accreted and 
ejected core material, which corresponds to an outflow velocity equal to the escape velocity of the collapsing core. For various 
values of the core flattening parameter ke, the fiducial value X — 0.32(1 + Hn) (equation circles) leads to 0csc — 67° from 
the axis. Bottom: The efficiency £coro of core collapse (equations |^ and solid lines) and an approximate formula (equation 
p^ , dotted lines). Curves are plotted for < kg < 10, corresponding to 1 < (1 + Ho) < 2, where approximation ( |33| ) is known to 
be accurate. 



Although the above argument assumes a constant 
value of Hq among cores, the result that the IMF 
is very close in slope to the core mass function is 
independent of this assumption. This follows from 
the fact that 0.3 < ecore < 1 (figure ||). The two 
slopes cannot differ greatly if the stellar mass is al- 
ways within a factor of three of the core mass. 

Features in the core mass function (such as a peak) 
appear in the IMF at a mass that is smaller by the 

factor ^core* 

5.4. Comparison with Observations 

We have adopted a scenario in which a protostel- 
lar core is destroyed simultaneously by accretion onto 
a central star (near the equator) and by protostellar 
outflow (near the poles). This is broadly consistent 
with observations that indicate both processes are 



important in dissipating cores, such as those of My- 
ers et al. (198^ ) and |Ladd, Fuller, k Deane (19980 . 



This geometry is confirmed by observations of the IRS 
2 outflow source in Barnard 5 ( fVelusamy &: Langei 
1998| ), which clearly discern an equatorial inflow co- 



existing with an axial outflow. 

Although we have not calculated the outflow struc- 
ture in detail, the current theory implies the existence 
of a massive, slow, conical outflow that reaches the 
core escape velocity at about 23° from the equator for 
the fiducial values of the parameters. This structure, 



which also appears in the theory of Li fc Shu (1996 ), 
is likely to correspond to the massive, poorly colli- 
mated, 1 km s~^ biconical outfiow observed around 



T Tauri by Momose et al. (1996| ). The outflow cone 
observed by Velusamy and Langer extends to ~ 27° 
from the equator, in rough agreement with our pre- 
diction. More speciflc comparisons would require a 
solution to the problem of simultaneous infall and 
outflow, which is beyond the scope of this paper. 

We derive values of Score that vary from 25% to 
75% for the fiducial value X = 0.32(1 + Ho), for val- 
ues of the core flattening parameter < kg < 10 (i.e., 
< Hq < 1). Thus, we expect that spherical cores 
will lose three-quarters of their mass while forming a 
star, whereas cores at the upper end of the acceptable 
range of flattening ( p^i fc Shu 1996| ) will lose only a 
quarter of their mass. The stellar initial mass func- 
tion therefore depends not only on the mass function 
for cores, but also on the relative degrees of flatten- 
ing of low and high-mass cores, and on changes of 
cores' line widths with mass. Observations that in- 
dicate a strong similarity between the core and star 
mass functi ons (in both shape and peak mass), such 
as those of Motte et al. (1998|) , require that Ecore be 



roughly constant and of order unity. This is certainly 
consistent with the above theory, since X{1 + Hq)^/"^ 
is typically of order unity and is small, so that ^core 
is indeed close to unity (eq. ^4|). Observations that 
indicate only similar slopes between the two distribu- 
tions (e.g., Testi fc Sargent 1998|) do not constrain the 
value of Ccore • They only require that it does not vary 
rapidly, a result that is consistent with the results of 
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6. INSTANTANEOUS STAR FORMATION EFFICIENCY 
IN CLUSTER-FORMING CLUMPS 

Clumps massive enough to form clusters of stars 
differ from the protostellar cores considered in the 
last section in several key respects. First, it is un- 
likely that there is a strong correlation between the 
geometry of the clump and the direction of the out- 
flow. Provided the clump is not highly flattened, we 
can obtain a reasonable approximation to the effect 
of the wind by assuming that the clump is spherical 
[Qifi) = !]• As a result, the star formation efficiency 
is given by equation (27). Second, much of the mass 
of the clump is generally not involved in the forma- 
tion of the star. An upper limit on the mass of a 
star that forms in the clump is given by the result 
for a protostellar core, in which all the mass is either 
accreted or ejected: 



for clumps, the efficiency parameter becomes 



*, max — £^corc-^cl 



2XMci 



1 + [1 + 4X2(1 + /^)2]l/2- 

There is no restriction on the formation of stars less 
massive than this, but stars with m^^ < rrii^fi have out- 
flows too weak to eject any matter from the clump. 

A further difference between clumps and cores is 
that since the clumps contain the cores, the cores 
are generally significantly denser than the clumps in 
which they are embedded, and their free-fall times 
are correspondingly shorter. Insofar as the accretion 
time (and therefore tw) is tied to the free-fall time, as 
in inside-out collapse models for star formation ( Shrj 
1977|]|, it follows that the star formation occurs on a 
time scale shorter than the free-fall time of the clump. 
In this case, the effect of gravity on the shell's motion, 
which is characterized by the factor Cg (eq. |A13 ), is 
relatively minor: 



9 - Sfep 
8 - 3kr. 



1/2 



(52) 



which is of order unity for kp < 2. We shall generally 
adopt kp = 1 for clumps, which is similar to or some- 
what shallower than the clump structures inferred 
from observations (for instance, Yun fc Clemens 
19910 estimate 1 < fep ^ 2 in Bok globules). 

The final distinction between clumps and cores is 
that the clumps are generally not isothermal, and 
the rms velocity is often significantly greater than 
the thermal sound speed. Since both cores and 
star-forming clumps are gravitationally bound, it fol- 
lows that the escape velocity is generally greater for 
clumps than for cores. Normalizing Vesc to 2 km s~^ 



X = 0.29 



ln(2/0i 



0) 



In 200 V 2 km s 



40 km s 



For X = 0.3 and fw = 0.2, we have m^, 



(53) 



, = 0.27Mei 

from equation (p|). Thus, for M^i ^ 500 Mq, a clump 
has enough mass to make a star with a mass up to 
n^*,ma.x ~ 140 Mq, significantly above the upper mass 
cutoff in most star clusters in the Galaxy, even those 
that are far more massive than the clumps we are con- 
sidering ( pVIassey et al. 1995| ). Thus, we do not expect 
the formation of a single massive star to consume all 
the gas in a star-forming clump. 

We shall neglect the loss of mass due to the winds 
themselves [m^(l — (pw) ^ M^j] for star formation 
within a clump. This approximation is justified be- 
low by the fact that M^j ^ for typical conditions 
(m^ ^ rri-i^fi). Furthermore, the accuracy of this ap- 
proximation is improved somewhat by the fact that 



(pw is typically about 0.5 in clumps (Matzner 1999). 
This implies that the error in e introduced by ignor- 
ing the contribution of niw to Mej is only a couple 
percent when e ^ 50%, as we will find is typically the 
case in regions without massive stars. 

To estimate e for a star- forming clump, we assume 
that all the star formation occurs at the center of the 
clump. As we shall see below, this apparently dras- 
tic approximation is in fact reasonably accurate. The 
critical angle for the escape of the swept-up shell is 



1 



nl/2 



(54) 



from equations (|T^) and ([25|), provided, of course, 
that > m^^Q. For m^^max ^ "rn* ^ ^*,o, as is 
typically the case, we have ^esc — 1 and equation 
(P^) for the star formation efficiency simplifies to: 



e + ■ H2/9o) [ v,sc / Uokms-i 



(55) 

Because we do not expect the intrinsic properties of a 
protostellar wind, and fwVw, to correlate with the 
escape velocity of the protostellar clump, equation 
(p^) indicates that the efficiency with which clusters 
of low-mass stars are formed is a simple, increasing 
function of Vesc- An escape velocity of order 2 km s~^ 
leads to e ~ 40%. 

Equation (|55| ) assumes that the mass ejected by 
each outfiow can be tallied independently of other 
outfiows. This can be verified by considering two 
outfiows that occur simultaneously, which maximizes 
their interaction. If the angle between their axes 
is significantly greater than ^esc for either, then the 
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momentum output from each is negligible within an 
angle 9esc of the other's axis. Then, Mgj is given 
by equation (^) for each outflow individually. Al- 
ternatively, consider the case in which the outflows 
are completely aligned. They would then act as a 
single outflow whose momentum (and wind mass, 
assuming identical values of f^) is the sum of its 
two components. So long as the combined mass of 
stars that forms simultaneously is still much less than 
max — Mci/3, so that our assumption 1 — /^esc ^ 1 
still holds, equation ( |26| ) shows that Mgj is then the 
sum of the values it would take for each outflow indi- 
vidually. In either case, and thus quite generally, the 
relative timing and orientation of outflows is of little 
importance when estimating the mass they eject. 

6.1. Numerical Models: Off-Center, 
lonization-Regulated Star Formation in 
Clusters 



We wish to assess the accuracy of equation ( pSD , in 
which we assumed that each star forms at the exact 
center of its clump. In this section, we present calcu- 
lations of £ for off-center star formation. As above, we 
assume that the clump is spherically symmetric with 
Pc\ oc r~^i'; we set kp = 1. We draw all stars from 
the same initial mass function independently of their 
position within the clump, and orient their outflow 
axes at random. The IMF is taken to be a Miller- 
Scalo log- normal distribution, with a lower cutoff at 
0.057 M0, so that = 0.5 M©. We take stars to 
be formed in the clump at a local rate t~} = M^jM^x 

equal to the ambipolar diffusion rate t~^, as in th e 
photoionization-regulation model of McKee (19891) . 
Because ^ad oc if cosmic-ray ionization prevails, 
we take the star formation rate to vary with the free- 
fall rate in the region of high extinction. The am- 
bipolar diffusion rate decreases sharply due to far- 
ultraviolet ionization as the visual extinction to the 
surface Ay^ falls below about four magnitudes. Al- 
though this transition can be described analytically 
( pVIatzner et al. 20001 ), for simplicity we shall take the 
transition to be sudden: tg^(r)~^ oc t^(r)^^ oc p{ry/'^ 
if Avs > 4 and t~^ = if Ays < 4. 

For each star, a grid of solid angles is constructed 
around its outflow axis. The wind momentum is 
distributed on this grid, whose resolution is concen- 
trated toward the axis for improved numerical accu- 
racy. The clump mass per solid angle is also com- 
puted, using the distance from the newborn star to 
the center of the clump, and using the angle between 
its axis and the radial direction. The ratio between 
dpw/dVt, and dMd/dQ gives the shell velocity at the 
clump surface; this is compared to CgVesc, and the to- 
tal ejected mass is tallied. Thus, the effect of gravity 
is included as an overall factor, rather than consid- 



ered in the motion of each outflow element. 

In Figure |2|, we compare the results of this nu- 
merical investigation to our analytical formulae. Al- 
lowing for distributed star formation and averaging 
over the initial mass function changes the estimated 
efficiencies appreciably only when the assumption 
m^^o m^^max breaks down (near the thin 

lines that represent these limits in Figure §). For 
instance, the curvature of the numerical contours of 
£ away from the analytical contours at the bottom 
of the plot can be attributed to the smallness of 
JTT-*, max/^* in that region. We have also compared to 
a calculation of e that averages over the stellar initial 
mass function but assumes each star forms at the cen- 
ter: this reproduces the off-center calculation almost 
exactly, except for m^^o > In general, assuming 
that each star forms at the center yields an excellent 
approximation of the star formation efficiency. 

6.2. Comparison with Observations 

The above theory implies that the clumps that give 
rise to clusters of low-mass stars, like those in L1630 
{nn ~ 10^ cm-3, 200 < < 500 M©; [Lada 1992|) , 
form each star at an efficiency 30% ^ e ^ 40%. This 
range of efficiencies agrees with the theoretical argu- 
ments by Mathieu (1983 ) and Elmegreen "(198^ ), who 
compared the properties of known embedded and re- 
vealed bound clusters of similar masses and deter- 
mined that star formation efficiencies of about 30% 
or 55% would be necessary for the embedded bound 
clusters to evolve into the revealed bound clusters, in 
the respective limits of slow or rapid gas dispersal. 

As described in the observable efficiency is not 
the instantaneous efficiency, e, but instead the ratio of 
stellar to total mass, SFE. Whereas e is a simple func- 
tion of the current escape velocity from the system 
(eq. [55| ), SFE is expected to vary from zero to unity 
as stars form and gas is dissipated; e is only a rough 
estimate for the characteristic value of SFE in the 
midst of this process. Without knowing the evolution 
of the system, it is therefore difficult to compare an 
observation of SFE to our prediction of e. In L1630, 
the observations of |Lada (1992]) imply SFE ~ 7% 
(NGC 2023 and LBS 23), 19% (NGC 2071), and 42% 
(NGC 2024 and 2068) if one assumes each detected 
source signifies 1 Mq of stellar material. This as- 
sumption is relatively uncertain due to variable ex- 
tinction in the area, luminosity variation in the em- 
bedded sources, and the completeness limit of the 
survey. If instead each source signified the mean stel- 
lar mass (~ 0.5 Mq), the inferred efficiencies would 
be about a factor of two lower. A direct comparison 
with theory is therefore difficult for both theoretical 
and experimental reasons, but it is encouraging that 
there is rough agreement for the richer clusters, which 
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Fig. 2. — Estimates of the instantaneous star formation efficiency e due to protostellar outflows. Dashed lines: ana- 
lytical estimate, equation (|55|). Solid lines: Numerical evaluation of the same quantity in a protostellar clump with the 
structure p oc and gas conversion rate scaling as p^^^ in regions extinguished by more than four magnitudes from the 
surface. In this calculation, the effect of gravity on outflow motions has been included with an overall correction factor 
Cg (equation |]5^ ). Dotted lines indicate the mean extinction of the clump. Above the dash-dot line, stars more massive 
than 20 Mq are likely to form; our estimates are therefore only rough upper limits in this region. Below the lower thin 
line, m^^ max is less than twice the mean stellar mass; above the upper thin line, m^.o exceeds the mean stellar mass. The 
analytical approximation requires to*,o 'tT* ^7i*.max, and breaks down near the thin lines. 



may be more evolved. 

(Jnishi et al. (1998) have conducted a survey in 
Taurus for dense gas and young stellar objects in var- 
ious phases of evolution, in order to ascertain the evo- 
lutionary sequence leading to cluster formation. They 
find that only those clumps (and all those clumps) 
with mean visual extinction Ay > 8 mag are actively 
forming stars, in agreement with the theoretical pre- 
diction of McKee (1989| ). However, some clumps of 
lower extinction contain more evolved objects that 
presumably formed while the column density was 
above this threshold. Considering the decrease in col- 
umn density to correspond to a loss of gas mass, and 
multiplying by the mean area per star, these authors 
estimate that about 4 Mq is lost per observed object, 
i.e., e(obs) ~ 11% if each object represents 0.5 Mq 
of stellar material. This observational determination 
of efficiency is also quite uncertain. It should be con- 
sidered a lower limit for several reasons: the stellar 
sample is likely to miss a significant fraction of weak- 
line T Tauri stars ( Bricefio et al. 199E , 199§| ); stars 
may have drifted away from the older clumps; and 
the decrease in column density could be partly due 
to an expansion of the clumps over time. Moreover, 
the currently star-forming clumps have not created 
their full complement of stars: on average, one might 
expect about half the stars to have been born by now, 
in which case e(obs) ~ 22%. The star-forming clumps 
in this sample have escape velocities ranging from 0.8 



to 1.4 km s^^, for which we estimate e ~ 25ib5% from 
equation (|55D. Again, the theory predicts efficiencies 
similar to the values inferred from observation. 

7. CONCLUSIONS 

The star formation efficiency is a fundamental 
property of the process of star formation. Star for- 
mation is often said to be inefficient because it occurs 
much more slowly than it potentially could, as dis- 
cussed in the Introduction. However, when defined 
in terms of the fraction of the available mass that 
goes into stars, it is rather high — at least for low- mass 
stars. For such stars, the efficiency is determined by 
the powerful protostellar winds that accompany the 
formation of the star. These winds produce a to- 
tal momentum per unit mass of star formed of order 
rUiuVw/^* = fw^w — 40 km s""*^. If this momentum 
were transferred to the surrounding gas isotropically, 
then the mass ejected would be much greater than the 
mass that goes into the star, and the efficiency would 
be low (Nakano et al 1995). In fact, however, pro- 
tostellar winds are highly collimated. To determine 
the star formation efficiency resulting from collimated 
outflows, we assumed: (1) the interaction between the 
protostellar wind and the ambient gas is mediated by 
thin, radiative shocks that lead to radial outflow; (2) 
except possibly near the wind axis, the mass of the 
wind is small compared with the mass swept up by 
the wind shock; and (3) there is an angle ^esc relative 
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to the wind axis inside of which all mass is ejected 
and outside of which mass can accrete toward the 
central source, with a fraction /^/ (1 + fw) going into 
the protostellar wind. Assumptions (1) and (2) are 
consistent with observations of the momentum dis- 
tribution in protostellar outflows ( Matzner fc McKeg| 
1999|]. 

We additionally assumed (4) that fwVw is the same 
among stars of different masses. This assumption, 
which is justified by the fact that deuterium burning 
regulates protostars' escape velocity QStahler 1988|) 
and thus the wind velocity ( ^hu et al. 1994 ), leads to 
an estimate of the star formation efficiency that is the 
same for most stars and depends only on the escape 
velocity of the region (eq. |5^). If this assumption 
did not hold, e would vary among stars. However the 
conclusions of § ^ would still be valid, as the star for- 
mation efficiency for a clump is an average over the 
IMF. 

The principal effect of wind collimation on the mo- 
mentum transfer to the ambient gas is to concentrate 
the momentum in a mass that is smaller than in the 
isotropic case by a factor of order l/ln(2/0o)i where 
$0 ~ 10~^ is the size of the central region of the out- 
flow (eq. |ll]). Gas is ejected if its velocity exceeds the 
escape velocity, Vesc', as a result, the mass ejected is 
given by Mejfesc ~ "i^/u>t'u>/ In (2/6*0). If the ejected 
mass is larger than the protostellar mass, then the 
star formation efficiency is approximately 



3ln(2/go) 

fw^w 



(56) 



This result, which applies both to the formation of 
individual stars in a core or the formation of a clus- 
ter of stars in a star- forming clump, shows that the 
efficiency of formation of low-mass stars is 

- independent of the mass of the star; and 

- only weakly dependent on the properties of the 
core or clump, through fgsc oc (Md/Rci)^^'^ oc 

Consequently, the efficiency is only weakly dependent 
on the geometry of the star-forming region and the 
location of the outflow source within it, as the re- 
sults of §§^ and 6.1 demonstrate. These conclusions 
are i n quali tative agreement with those of Nakanc| 
et al l (19951) , but our inclusion of the collimation of 
the outflows increases the estimated star formation 
efficiency by about an order of magnitude, so that 
e ~ 0.25 - 0.75. 

We have evaluated the star formation efficiency 
separately for protostellar cores and for star-forming 
clumps. Protostellar cores are expected to be flat- 
tened in a direction perpendicular to the outflow; 



since the outflow and inflow are in different directions, 
the star formation efficiency is increased over the case 
of an isotropic ambient medium. We evaluated this 
effect quantitatively using the Li-Shu (1996) model 
of isothermal, magnetized protostellar cores (see Fig. 
1). Because e is independent of the stellar mass and 
only weakly dependent on the core mass, we predict 
(in agreement with Nakano et al 1995) that the initial 
mass function is just slightly flatter than the mass 
distribution of protostellar cores, with a power-law 
index (0.9 — l)acore- 

Our model differs from most previous proposals in 
that the wind acts to limit the fraction of the mass 
that can accrete, but does not cause accretion to shut 
off entirely. This is consistent with the fact that ac- 
cretion is necessary to drive a powerful wind. Instead, 
accretion terminates because only a finite reservoir 
of gas is available. In contrast to the proposal of 
Velusamy fc Langer (1998| ), for instance, the angle 
dividing accreted and ejected material takes a charac- 
teristic value during the main accretion phase. Varia- 
tions of ^esc may result from a decline in the accretion 
rate, but do not cause it. 

Another distinction of the theory presented here is 
that the division between accretion and ejection oc- 
curs in swept-up material traveling outward at the 
escape velocity - material outside the region of sig- 
nificant infall and far outside the centrifugal radius. 
Therefore the complicated interaction among outflow, 
infall and rotation, which determines the shape of the 
outflow cavity at its base, does not affect the effi- 
ciency. All of the material at the base of the flow is 
destined to accrete, for rotational flattening and infall 
only act to lessen the effect of the wind. 

We have found that the star formation efficiency 
is controlled by the efficiency parameter X (equa- 
tion Q), and in the case of magnetically flattened 
protostellar cores, also by the flattening parameter 
(I + Hq). The factors that introduce uncertainty into 
the theory are primarily those that affect X. For in- 
stance, the factor Cg that corrects for the effect of 
gravity on the outflow motions is well-understood for 
clumps of hydrogen number density below 10^ cm~^, 
in which winds are impulsive. As we discuss in 
we expect that winds are approximately impulsive for 
higher mean densities as well; however, the accuracy 
of this approximation is difficult to assess. For the 
formation of an individual star, the wind duration is 
longer than the free-fall time of the core. This raises 



the question, addressed in |A.2| of the Appendix, of 
what happens to material when it reaches the core's 
surface. The results presented in ^presume that the 
shocked shell decouples from the wind outside the 
core; however, equation ( [A17| ) shows that Cg would 
be three times lower in the opposite limit of perfect 
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coupling. 

A second source of uncertainty is the parameter 
fwVw that relates a star's final mass to the momen- 
tum of its wind. As discussed in our value, 
fwVw = 40 km s~^, is consistent with current obser- 
vational and theoretical estimates; however, both are 
somewhat uncertain. 

Our estimate of the effective broadening angle, 6q, 



of protostellar winds is also not precise. Matzner & 
McK ^e (1999| ) estimate — 10 ^ by comparing the 
mass-velocity distributions of outflows to the predic- 
tions of a simple model for their dynamics, and by 
requiring that the model fit observations for realis- 
tic inclinations of the wind axis to the line of sight. 
Note that enters logarithmically in the efficiency 
parameter X, as log(2/0o); therefore, we consider it 
a minor source of uncertainty as compared to fwVw 
However, because the intensity of protostellar winds 
along their axes is proportional to 6q'^, this parameter 
plays a more important role in determining the min- 
imum stellar mass q capable of ejecting any mass. 
The efficiency increases dramatically (and equation 
[SStl breaks down) for clump masses large enough that 
nii^fi becomes comparable to the mean stellar mass. 
For 6*0 = 10~2, this happens when Md ;^ 10"^ M© 
(the upper thin line in figure ^). Another condition 
for the validity of our theory is that massive stars do 
not dominate the ejection of mass; as shown by the 
dash-dot line in that figure, we expect that to happen 
when Mci ^ lO^-^ M©. Thus for 6*0 = IQ-^, the for- 
mer is the more restrictive criterion; if were greater 



than about 10~^'^, the latter would be. 

Lastly, we expect that the assumption of radial mo- 
tion is not perfect. Most likely, this implies that our 
estimate of X should be adjusted by a correction fac- 
tor of order unity. 

For star-forming clumps, we found that the star- 
formation efficiency is typically somewhat less than 
0.5, although slight variations in the parameters can 
bring e above 0.5. As discussed in §^]^, such high ef- 
ficiencies might be expected to leave the cluster grav- 
itationally bound. To show why most clusters are 
unbound, we must consider the dynamics of cluster 
formation, a topic we address in a future paper. 
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APPENDIX 

EFFECT OF GRAVITY ON OUTFLOW MOTIONS 

Because fesc marks the boundary between ejected and retained material, and because gravity is important for 
the dynamics of a shell if its mean velocity is about Vesc, we wish to refine our estimate of the escape criterion 
accounting for gravity. 

Assume that a star forms at the center of the clump, and that the clump's gravitational potential is not 
perturbed by the passage of the outflow. Take the distribution of momentum injection to vary with angle away 
from the outflow axis as 

dpnj _ PwPjl^) .... 

In - 

our model for P(/i) is given in equation (p!o|). We take the distribution of clump density to be 

Pel = PoiQ{p)r-^' (A2) 

for < r < i2ci) where the normalization of the angular distribution of the density is /q^ Q(/i)ci// = 1. The mass 
per unit solid angle enclosed within radius r is therefore 

= ^Md(r)g(/z) = "-^Q{^^y-'% (A3) 

where the total mass inside r is Mc\{r) = CMr^~^'', and cm = 47r/Joi/(3 — kp). 

To determine the condition for the wind to eject material from the surrounding clump or core, we make 
several approximations (§01): We assume radial, momentum-conserving flow, and we neglect dniw/dVL compared 
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with dM(.\/(Kl in regions in which the velocity of the sheh is small compared with the wind velocity, Vs <^ Vy,. 
In addition, we adopt the monopole approximation for gravity, which is generahy quite good, even for flattened 
distributions ( Li fc Shu 19971 ). The equation of motion for a sector of the sheh is then 



d dMci 
lit dn 



dpw 

dn 



GM dM 



or 



r2 dn ' 



(A4) 



(A5) 



where we have divided through by Q{^)- In the monopole approximation, the kinematics of the shocked sheh 

depend only on the ratio P{fi)/Q{fi), not on either one independently. 

If we neglect the effect of gravity on the propagation of the shock, as in the main text, equation ( A5) integrates 
to give 

^^^\pw{t). (A6) 



MclVs 



Usually, we shall adopt the criterion that a sector of the shell will escape from the clump if its velocity exceeds the 
escape velocity. The critical angle at which escape is marginally possible, ^esc = cos~^ ^esc, is then determined 
by setting Pw{t) = pw, the total wind momentum, and Vg = fesc in equation ([A^). Including the effect of gravity 
on the dynamics of the shell increases the required momentum by a factor Cg >\ that we shall determine below: 



zPw = CgMclVe 



(A7) 



To solve for the dynamics of the shock in the clump including the effects of gravity, we express Mci(r) as 



CMr' 



3 kp 



and d/dt as Vgd/dr. Equation ( |A4| ) then becomes 



d 

CMVs-rVsr 
dr 



3 kp 



(A8) 



We now introduce the variables x = r^ and y = 



,3-fc 



in terms of these variables, equation (A8) becomes 



- kp dy^ 



:Pv 



Integrating, and transforming back to radius and velocity from x and y, we find 



P{^^) 



PwT 



MdVescirf , C 



+ 



Mel' 



where we have defined Vesc{r) = ^j2GMc\{r) /r as the local escape velocity [so that fcsc(^ci) = 'I'c 
constant of integration. 



(A9) 

(AlO) 
and C is a 



Impulsive Winds 

Consider the case that the wind momentum is injected impulsively, i.e, tyj ^ %. Then, Pw = ^ during the 
shell's expansion, so that (multiplying through by Mf.\) 

M>esc(r-)2 



Mlvl 



8 - 3A;„ 



C. 



(All) 



Evaluating this equation at r = 0, where Mcifesc('^) = 0, we find MciQ{lJ,)vs = PwPilJ-)', therefore C 
[p^P(^)/Q(m)]2, and 



PM. 



8 - 3ko 



The condition for the escape is ^^(-Rci) > Vesc, or 



PM 



Pw > 



9-3kp 
8-3kr. 



1/2 



Mc\Ve 



(A12) 



(A13) 



Thus, in this case the factor by which the critical wind force must be enhanced to offset the effects of gravity 
on the dynamics of the shell is Cg = [(9 — 3kp)/{8 — 3kp)]^^'^ (cf eq. A7). 
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Steady Winds 

Next, consider what happens if the wind momentum injection is constant throughout the expansion of the 
shell, so that pw maintains the steady value pw/tw Furthermore, let us suppose that the wind continues to 
inject momentum after the shell breaks out of the clump. Then, a naive estimate for the criterion for escape is 
that the wind force overcome gravity at the clump surface, 



— — Mcl^esc- 



(A14) 



This estimate gives Cg = (7r/4)(t^/tff ) for the factor by which the momentum in the wind must be enhanced 
to overcome the effects of gravity on the shell dynamics. This factor is large when tw ^ since only the 
momentum delivered in a time of order is effective in ejecting the shell. In fact, the condition Cg > 1 sets a 
lower bound on the value of tw/tff for a steady wind. 

Consider the dynamics of the outflow shell after it has broken free of the clump, so that its mass is constant. 
Then, since we are assuming that p^ is constant in time, we can integrate equation ( [A5D directly, using 
d/dt = Vgd/dr: 

E{r) - E{R,,) = " ^ci), (A15) 

where E{r) = Mc\v^/2 - GM'^Jr is the energy of the shell and E[Rc\) is its energy at the clump surface. If the 
wind were no longer blowing, the shell's energy would be constant and positive energy {vl > GMl/r) would 
lead to escape. However, we are assuming that the wind continues to blow well past the time that the shocked 
shell reaches the clump surface. In this case, the shell might escape even if it has negative energy, because of 
the wind force. (Below, we shall consider the case that the outflow decouples from the wind after breaking 
out.) 

We shall evaluate Vs{Rc\) using equation ( AlO] ), and then use equation ( A15| ) to determine the criterion for 
escape. We readily see that C = for a steady wind, for all the other terms in equation ( AlO ) are zero at the 
origin; therefore. 



2 p^rPiii) M^iVescirf 



4- kp twQ{n) 8 - 'ikp 



(r < R, 



■clj 



Using the value of E(Rci) from this in equation (|A15| ), and requiring that dE{r) 
Vg{r) > 0] for all r, we find the critical value of the wind force: 



> 



Pw > 



1/2- 



X ——MclVe 



(A16) 

-GMcidMci/r [so that 
(A17) 



This more precise wind force criterion is 35% and 31% lower than the simple estimate of equation (A14) for 



kp = and hp = 1, respectively. The two agree exactly for kp = 2, because gravity and wind force scale in the 
same manner inside the clump in that case. This implies that the shell velocity at the surface vanishes if the 



wind force equals the critical value if fcp = 2, exactly as we assumed in estimate ( |A14 ) 



Equations ( |AT4D - ( [ATtD assume that the wind continues to deposit momentum into the shell efficiently, even 
after the outflow has emerged from the clump. There are two reasons that this assumption may be incorrect, 
which make the above estimates lower limits for the critical wind force. First, if is not too much longer 
than iff, there is the possibility that the wind force will taper off before the shell has achieved escape velocity. 
Second, there is the possibility that the shell will be disrupted by a thin-shell instability. The former effect 
can be ignored, because it is only important when — tg. The latter is more serious, as it would remove the 
supporting effect of the wind (partially or completely) once the instability set in. 

Shock-bounded thin shells are subject to nonlinear instabilities (e.g., Vishniac 1994 ); these are expected to 
become stronger if the dense fluid (the ambient gas, in this case) is "upward", i.e., if the shell accelerates, or 
at least decelerates more slowly than gravity would predict (E. T. Vishniac, private communication). Writing 
out d{dMci/dQ) /dt as r'^VsPd and dp^/dQ as /j, 
the actual acceleration and gravity is 



we find from equation (|A4| ) that the difference between 



dt r2 



dMci/dn 



Pc\Vs 



(A18) 
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Since gravity is an additive term in the acceleration, we see that the shell decelerates more slowly than gravity 
would predict (and is unstable) if the same shell would accelerate, were gravity suddenly turned off. Equation 
( |A18 ) shows that this only occurs if the ram pressure in the wind exceeds the ram pressure of the ambient 
gas. For a steady wind, this requires that the ambient medium must steepen so that kp > 2 (declining winds 
would require even steeper profiles). We therefore expect that shells might become partially decoupled from 
their winds as they emerge from a clump, not before. 

Although this instability is nonlinear, and may be suppressed by magnetic fields in the outflow shell or in the 
wind, we can obtain an upper limit to the critical wind force if we assume that outflows decouple entirely from 
their winds after breakout. Then, the criterion for ejection is once again Vs{Rc\) > v^sc- From equation ( A16| ), 
this requires 

T^rr-rPw > (4 - /c„) ^ x ^ = (4 - k^) x — Md^esc- (A19) 

The factor Cg in this case is the coefficient of MdUesc in the equation above. We see that the wind force required 
to eject the shell is much stronger if the shell decouples from the wind after breakout than if the wind continues 
to provide support, the factor between equations ( A17 ) and ( |A19 ) being seven or three if kp is zero or two, 
respectively. 
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